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We describe the creation of nonclassical states of microwave radiation via ideal dichotomic single 
photon detection, i.e., a detector that only indicates presence or absence of photons. Ideally, such a 
detector has a back action in the form of the subtraction operator. Using the non-linearity of this 
back action, it is possible to create nonclassical states of microwave radiation, including squeezed 
and cat-like states, starting from a coherent state. We discuss the applicability of this protocol to 
current experimental designs of Josephson Photomultipliers (JPMs). 



The generation of nonclassical states of light is an 
important test of the foundations of quantum mechan- 
ics and a necessary precursor to implementing quantum 
communication and computation protocols in many ar- 
chitectures [1-3]. While the methodology for creating 
nonclassical light at optical wavelengths has been stud- 
ied extensively [4-6], the technology to create quantum 
states with larger and larger wavelengths has recently 
become available with advances in circuit-QED. 

In this paper we present a novel way to generate non- 
classical states of microwave radiation in a long wave- 
length transmission line using only detection by an ideal 
binary d(rtcctor. the Josephson Photomultiplicr (JPM). 
The protocol only involves radiating a microwave cavity 
with coherent light, and post selection based on single 
photon detection by JPMs without further manipulation. 
In addition, our protocol applies to any detection mech- 
anism with a back action resembling that of the subtrac- 
tion operator (equation (1)) and so can be generalized 
to other quantum systems, in particular, other supercon- 
ducting circuits where strong photon-detector coupling is 
possible [7-9[, but also at optical frequencies. 

In the microwave regime, squeezed states [10-12] and 
cat-like states [13[ of microwave radiation have been gen- 
erated through the QND measurement of a microwave 
cavity with a Josephson parametric amplifier. We con- 
sider a projective detection scheme whereby the isolated 
cavity is left in a pure state, a protocol more similar to 
the Quantum Zeno detection scheme proposed for cavity- 
QED experiments [14]. We also note that the creation of 
arbitrary nonclassical states via the controlled interaction 
of a superconducting qubit with a microwave cavity has 
been achieved experimentally ]8, 9]; here we show that 
a certain class of nonclassical states can be constructed 
more efficiently than thought previously. 

The JPM, a current biased Josephson junction related 
to the phase qubit, has been shown experimentally [15] 
and theoretically [16-18] to be an effective single mi- 
crowave photon counter. Previously, we have shown that 
for a JPM under optimal conditions the back action of 
photon detection is the photon subtraction operator [16], 
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a non-linear operator that can be related to the photon 



lowering operator by d = B\fti, but cannot be expressed 
as a linear combination of photon creation and annihila- 
tion operators. Also, note that B is not invertible, and 
hence not unitary. 

In this paper we show how to use the noncommunta- 
tivity of the detection back action with coherent displace- 
ment pulses to achieve single mode quadrature squeezing 
of an input coherent state as well as other nonclassical 
states, including generalized Schrodinger cat states [4]. 
We also discuss the success probability of our protocol 
under realistic operating conditions. 

The generation of squeezed states of microwave ra- 
diation using JPMs follows a simple protocol. The 
cavity is initially prepared in a coherent state, |a) = 

J2^=o 7^1"') where a = jaje*'''", and is coupled 

to one or more detectors, each acting with back action B 

on the cavity after a photon is detected. After N photons 
are counted, a displacement pulse is applied such that the 
state is centred around —a in phase space. After N fur- 
ther photon detections are observed, the resulting state 
is a squeezed state. The optimal choice of N will be 
discussed shortly. 

Starting from the coherent state input, the probability 
for N photons to be detected is 

^^-^ ^ -1 r(iv) 
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where r(A'^, is the upper incomplete gamma function 
of N and \af [19]. ~ |ap^/iV! as |a| -S> 0; however, 
Pjv increases monotonically with |q!|, so the success prob- 
ability can be made arbitrarily close to unity with high 
power coherent pulses. 

It is straightforward to calculate the normalized post 
measurement cavity state after N detections, 

(3) 

and the average photon number n' = (a^a)^, is given by 

n = ^ ^ ^ (4) 
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or calculated straightforwardly numerically. After N de- 
tections, the next step is to displace the state by an 
amount a' ~ —y/n'e^'^" — a, so that the resulting state 
will be centred in phase space around —a. For this dis- 
placed state, N photon detection events will occur with 
probability [27] , and the renormalized cavity state 
will have the form 
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We will now show that the state p" is a squeezed state. 

To quantify the amount of squeezing, we calculate the 
variance of the squeezed quadrature 



Ap2 = TY[p"f] - Tr[p"j5]2. 
The quadrature observable p{4ia) is defined by p{ 
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where a is the annihilation opera- 



tor for the cavity microwave mode. The phase shift e~"^° 
accounts for the fact that this protocol squeezes along the 
phase space axes defined by the phase of the input co- 
herent state. Anything less than Ap^ = ^ indicates a 
squeezed state. The amount of squeezing is expressed in 
dB, by calculating 



S'(Ap) = lOlog 



10 



Ap2 



Ap; 



,2 

norm 



lOlogio (2Ap2 



(7) 



In addition, we can calculate how far the state p" de- 
viates form a minimal uncertainty state by calculating 
AxAp (where ^ = ^ (ae^*'''° -I- a^e*'''") is the conjugate 
observable to p). Figures 1(a) and 1(b) show 5'(Ap) and 
Aa;Ap respectively as functions of a and the number of 
detection events on either side of the displacement, TV. 

As can be see in figure 1(a), the maximum amount of 
squeezing possible on a given input state \a) increases 



monotonically with 



proportional to the power of 



the input pulse. Interestingly, for a given a there exists 
a finite N that achieves a global minimum in Ap^. One 
would be tempted to use this value of N in the protocol 
to create p" , however, as can be seen in figure 1(b) there 
are other concerns. 

As figure 1(b) shows, for a given a, AccAp is not mono- 
tonic in N . Since we want p" to be as close to a minimal 
uncertainty state as possible, while still maintaining a 
significantly squeezed quadrature, the optimal choice of 
N for the protocol would be at the AxAp local min- 
ima shown in figures 1(a) and 1(b) by the white curve. 
While this does not minimize Ap^ (and therefore maxi- 
mize squeezing) , it achieves a significantly squeezed state 
p" that is as close to being a minimal uncertainty state 
as is possible, which is what we consider optimal. 

Other nonclassical states of microwave radiation hav- 
ing = 211 /k rotational symmetry in phase space re- 
sult from generalizing this procedure (see figure 2). The 
squeezed state protocol discussed previously involves de- 
tecting the cavity state while centred at two points on 
a line through the origin of phase space, i.e., the case 
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FIG. 1: (a) S{Ap) and (b) AxAp of the state p" , as functions 
of a and the number of detection events, A^. ipa = for both 
figures. The white curve indicates the value of A'^ that gives a 
local minimum in AxAp, while maintaining a squeezed Ap^. 



k = 2. To generalize this protocol, we detect N photons 
at k positions equally spaced around a circle of radius \a\ 
in phase space, and take the first position on the positive 
real axis for simplicity. There are k steps to this gener- 
alized protocol, and for j = 0,1, ..k — 1, the jth step is 
to 1) displace so that the cavity state is centred around 
lale*^'^*, and 2) detect N photons. Finally, after detect- 
ing N photons at all k positions, we displace the state so 
that it's centred around the origin for comparison with 
the squeezed vacuum, and we refer to this final state as 
Pk- The detection stage at each position transforms pji 
to pj2 according to 



Pj2 



pj 



(8) 



where is the probability of detecting N photons from 
the state pji. The entire protocol will complete with 
success probability 



fc-i 



Prob(success) — JJ^ P^, 



(9) 
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which monotonically increases with \a\. 

As mentioned previously, the case k = 2 (ip^ = tt) 
corresponds to the creation of the vacuum with squeezed 
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quadratures. We find high (here and henceforth meaning 
above 99%) overlap with 

|*2> = Siz) |0) = e-K^(a*)'-^*(a)') |o) (10) 

by numerically searching over N and z, where S{z) is 
the ordinary squeezing operator with complex squeezing 
parameter z. In general, for k > 2, the states created by 
this protocol are very close to the analytic states 

l^fc) = S'^''\z) |0) = e-K^(a*)'-^*(")') |0) , (11) 

where S'^'^' (z) is called a generalized squeezing operator 
with complex parameter z [20, 21]. The first three oper- 
ators of this class are 
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S'-^\z)^S{z) 



(12) 
(13) 
(14) 



where S{z) is the squeezing operator of equation (10). 

Consider the final state after displacing to the origin 
Pk- To find z, we impose the condition that pf^ and \'i'k) 
have the same average photon number by setting 



(15) 



This results in excellent fidelity between the states, quan- 
tified by 



F[pfe]=Tr[pfe|*fe)(** 



(16) 



In fact, for |a| = 3 and A; = 2, 3, 4 the fidelity is > 99% 
in each case. The Husimi Q functions of these states are 
shown in figures 2 (a), (c) and (e). 

To verify that these states are indeed nonclassical, a 
suitable nonclassicality witness can be used. For this pur- 
pose, we use the entanglement potential of [22], defined 

by 



EP[p]^log2||aJ-||i, 



(17) 



where ap = [/bs (p |0) (0|) f/gg, for Ubs the unitary 
transformation of a 50:50 beam splitter, and aj^ is the 
partial transpose of ap [23[. A nonzero entanglement 
potential indicates that the state p is nonclassical, and 
indeed, as can be calculated numerically, the states cre- 
ated by the generalized protocol all have nonzero entan- 
glement potential (see supplementary material). 

Since the displacement stage of each step ensures the 
average photon number of each state pj is on the order 
of jap, each can be quite large, and as a result, the 
generalized protocol will have a significant probability of 
success. For example, for the modest \a\ = 3 of the states 
shown in figure 2, the success probabilities are all « 85%. 
Furthermore, the success probability will scale monoton- 
ically with |q;|, and so can be increased by increasing the 
initial input state power. 
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FIG. 2: This figure shows the Husimi Q functions for the 
states created by the generalized protocol, figures (a), (c) and 
(e), and those created by one additonal photon detection after 
displacement to the origin, figures (b), (d) and (f). 



This protocol can be generalized even further to create 
additional nonclassical states, including the squeezed odd 
Schrodinger cat state, with only a single additional detec- 
tion required. This extension to the generalized protocol 
is similar to that proposed in [6[ to prepare Schrodinger 
cat-like states, but does not require a beam splitter or 
photon number resolving detectors. 

If we remove one more photon from the final state pk, 
the resulting state has the form 



Pk 



BpkB 

^1 



(18) 



where is the probability of a single photon being de- 
tected from state pk- This procedure produces states of 
high overlap with 
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(19) 
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where now both z and /3 must be found numerically. 
These states all have non-zero entanglement potential, 
and are therefore nonclassical states of microwave radia- 
tion (see supplementary material). 

For k — 2, this additional detection will create a state 
very close to a squeezed odd Schrodinger cat state 



|vl/^>=5(2)(z) (1/3)- 1-/3)) 



For fc = 3 and /c = 4, we have 



\^',)^S^'Hz) (|/3) 
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(21) 



1*^) = 5W(z) (1/3) + |z/3) + |-;3) + . (22) 

The /c = 3 state is the state formed when the operator 
S'^'^\z) is applied to a Voodoo cat state [9]. The fc = 4 
state is the operator S'^^^z) applied to a coherent su- 
perposition of four out of phase coherent states, known 
as a compass state, which is known to have favourable 
decoherence properties [13, 24]. 

The Husimi Q functions of these squeezed cat-like 
states (for \a\ — 3) are shown in figure 2, all of which 
have greater than 99% fidehty with (20), (21), and (22). 
For the k — A case, it was necessary to slightly modify 
the shifts ai and 03 used in the generalized protocol by 
a small increment — |5|e*(^'^''+'''°) to achieve high fidelity. 
For the \a\ — 3 case shown, |^| « 0.1. 

The probability of successful generation of these 
squeezed nonclassical states is 



Prob (success)' = J| Pj 



(23) 



Unfortunately, these Pi are often very small, and for the 
squeezed cat-like states shown in figures 2 (b), (d) and 
(f) this results in a success probability less than 10%. 
However, as Pi grows monotonically with \a\, the success 
probability improves with higher power coherent pulses. 

In regards to experimental implementation of this pro- 
tocol with JPMs, superconducting microwave resonators 
are currently fabricated with Q-factors approaching 10^, 
which in the microwave regime will lead to cavity life- 
times on the order of lO^ns [25]. Thus, when JPMs with 
short T2 are used [16[, we can conservatively expect that 
as many as 10^ to 10^ measurements can be performed in 
the lifetime of the cavity. Practically, space requirements 
on a chip require a fast reset strategy for the JPMs, which 
is currently being developed [26[. 



In conclusion, we have shown how a combination of 
strongly coupled photon counting and coherent displace- 
ment can be used to create nonclassical states of radiation 
with high probability. This protocol can be realized in 
circuit QED using Josephson photomultipliers. 
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Strauch and R. McDermott. Research supported by 
DARPA within the QuEST program and NSERC 
Discovery grants. LCGG was supported by the Ontario 
Graduate Scholarship Program. 
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FIG. 3: This figure shows entanglement potential (EP) of the 
states created by the generalized protocol (for fc = 2, 3, 4) as 
a function of the number of detections, A'^, performed at each 
step in the protocol. The generalized squeezed states are the 
solid curves, and the dotted curves correspond to the cat-like 
states. These states have lal =3. 



The entanglement potential is not only a nonclassical- 
ity witness but is in fact also a nonclassicality measure. 
Here we examine how the nonclassicality of the states 
changes with the number of detections at each step of 
the generalized protocol. Figure 3 shows the entangle- 
ment potential of the generalized squeezed states and 
the cat-like states (for fc = 2, 3, 4) as a function of the 
number of detections at each step. As can be seen, for 
the squeezed state (k = 2), the maximum value of EP 
occurs around the optimal value of N determined in the 
main text. In addition, for each state shown in figure 3, 
the value of A'^ which gave the highest fidelity with the 
ideal state corresponds to a near maximum value of the 
entanglement potential. 
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